o 

(N 
>>■ 



Enhanced Group Analysis and Conservation Laws 
of Variable Coefficient Reaction— Diffusion Equations 
with Power Nonlinearities 

O. O. Vaneeva)^, A. G. Johnpillai^'^ , R. O. Popovych)^ and C. Sophocleous^^ 

^Institute of Mathematics of National Academy of Sciences of Ukraine, 
3 Tereshchenkivska Str., Kyiv-4, 01601 Ukraine 

Department of Mathematics, Eastern University, Chenkalady, Sri Lanka 
^ Department of Mathematics and Statistics, University of Cyprus, Nicosia CY 1678, Cyprus 
^ vaneeva@imath.kiev.ua, andrewgratienj@yahoo.com, ^ rop@imath.kiev.ua, ^ christod@ucy.ac.cy 

I A class of variable coefficient (l+l)-dimensional nonlinear reaction-diffusion equations of 

. the general form f{x)ut = {g{x)u"ux)x + h{x)u"^ is investigated. Different kinds of equiv- 

alence groups are constructed including ones with transformations which are nonlocal with 
respect to arbitrary elements. For the class under consideration the complete group classifi- 
cation is performed with respect to convenient equivalence groups {generalized extended and 
conditional ones) and with respect to the set of all local transformations. Usage of different 
equivalences and coefficient gauges plays the major role for simple and clear formulation of 
the final results. The corresponding set of admissible transformations is described exhaus- 
' tively. Then, using the most direct method, we classify local conservation laws. Some exact 

solutions are constructed by the classical Lie method. 



^ ■ 1 Introduction 

' In this paper we study nonlinear partial differential equations (PDEs) of the general class 

o 

:b 

with Lie symmetry point of view. Here / = /(x), g = g{x) and h = h{x) are arbitrary 
■ smooth functions of their variables, f{x)g{x) ^ 0, n and m are arbitrary constants. The linear 

d . case is excluded from consideration as well-investigated. We additionally assume the diffusion 

coefficient to be nonlinear, i.e. n 7^ 0. The case n = is quite singular and will be investigated 
separately. 

- Various simple transport models of electron temperature in a confined plasma are reducible 

r> I to equations of the form [19]. Here u represents the temperature, f{x) the density and 

' g{x) the density-dependent part of the thermal diffusion. The term h{x)u™' corresponds to the 

heat source. Furthermore a number of special cases of the class Q have been used to model 
successfully problems in mathematical physics, chemistry and biology [8,27,28,34,42]. 

While there is no existing general theory for solving nonlinear partial differential equations 
(PDEs), many special cases have yielded to appropriate changes of variables. Point transfor- 
mations are the ones which are mostly used. These are transformations in the space of the 
dependent and the independent variables of a PDE. Probably the most useful point transforma- 
tions of PDEs are those which form a continuous Lie group of transformations, which leave the 
equation invariant. Symmetries of this PDE are then revealed, perhaps enabling new solutions 
to be found directly or via similarity reductions. The classical method of finding Lie symmetries 
is first to find infinitesimal transformations, with the benefit of linearization, and then to extend 
these to groups of finite transformations. This method is easy to apply and well established in 
the last few years [5, 12, 15,30]. A large number of diffusion- type equations has been studied 
using Lie group analysis. See for example, [10,16,38]. 



f{x)ut = ig{x)u^u,), + h{x)u^ (1) 
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A goal in the present work is to classify the Lie symmetries of the general class (pQ). The 
problem of classification of Lie symmetries is dated back in the late fifties, where Ovsiannikov [31] 
determined all the forms of the well-known nonlinear diffusion equation 

Ut = [f{u)Ux]x 

that admit such symmetries. Since this latter work a number of articles on the classification 
of Lie symmetries of diffusion type equations appeared in the literature. See for example, 
[3, 7, 9, 10, 16, 25, 29, 38]. The Lie symmetries of ^ when f = g = h = 1 can be found in [9, 10] 
(see also [14, Chapter 10]). In fact, Dorodnitsyn [9, 10] carried out the classification for the 
nonlinear diffusion equation with a source. 



Recently Popovych and Ivanova [38] presented the classification of Lie symmetries for the general 
class of diffusion-convection equations 

f{x)ut = [g{x)A{u)ux]x + B{u)ux 

and later Ivanova and Sophocleous [16] derived the Lie symmetries for the class 



Furthermore in the present work, using the most direct method [39], we carry out two classi- 
fications of local conservation laws up to equivalence relations. Conservations laws of when 
f = g = h = 1 can be found in [11]. 

The presence of three functions in equation (Q) makes our task very difficult. However 
knowledge of the equivalence transformations for class enable us to consider an equivalent, 
but simpler form. For this reason, in the next section we construct different kinds of equiva- 
lence groups for the class under consideration. Extensive investigation and usage of equivalence 
transformations have the two- fold effect for further group analysis of class : reduction of the 
number of cases to be studied and simplification of form of both equations and basis elements of 
their Lie algebras in these cases. Moreover, presentation of classification results in a closed form 
becomes possible only after application of generalizations of usual point equivalence group. 

The group classification of the equations from class (0) is carried out in Section |21 It is 
shown that usage of the generalized extended equivalence group and conditional equivalence 
groups plays a crucial role in solving this problem. Importance of true choice of gauge for 
arbitrary elements with respect to equivalence transformations is discussed and illustrated by 
examples. Let us note that similar investigation on generalizations of equivalence group and 
different possibilities in choice of gauges was first fulfilled in [17,18] for class Q. 

In Section |1] we look for the additional equivalence transformations between obtained cases of 
symmetry extension, which are inequivalent with respect to transformations from the adduced 
equivalence groups. As a result, the problem of group classification in class with respect to the 
set of all local transformations is solved in passing. Connection between additional equivalence 
transformations and conditional equivalence groups is demonstrated. 

Results of Section 13 are key for this paper. The set of form-preserving transformations [22] 
(called also admissible [37] ones) of class is described exhaustively. Roughly speaking, any 
point transformation which links a pair of equations from the class under consideration is called 
an admissible transformation for this class. See [37] for rigorous definitions and problem state- 
ments. At the best of our knowledge, first the set of admissible transformations was described 
in non-trivial case, namely for a class of generalized Burgers equations, by Kingston and Sopho- 
cleous in [21]. Although the same problem has been already solved for a number of different 
classes of PDEs [6,22,23,37,40], these results are not well-known. 



Ut = [k{u)ux]x + q{u). 





(3) 
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The set of admissible transformations of class (pQ) has an interesting structure. The class (pQ) 
can be presented as a union of normalized subclasses. Some of the subclasses admit non-trivial 
extensions of the equivalence group, which will be conditional equivalence groups for the whole 
class (^. Moreover, there are subclasses intersecting each other. In such cases interaction of 
the corresponding conditional equivalence groups gives admissible transformations which can 
not be interpreted in group terms. By-products of investigation of admissible transformations 
are, in particular, exhaustive description of the additional equivalence transformations between 
symmetry extension cases and derivation of all condition resulting in non-trivial conditional 
equivalence groups. 

After presenting necessary notions and tools (Section inj, in Section d we classify local con- 
servations laws of the equations from class 

Studying equations (pQ) in the framework of Lie group analysis is completed in Section |H1 with 
construction of exact solutions for some equations from this class. We use results of fulfilled 
group classification and apply the classical Lie method. The main idea of the section is to show 
simplification introduced in finding exact solutions with equivalence transformations. 



2 Equivalence transformations and choice 
of investigated class 

In this section we derive the equivalence transformations of class (^. These transformations en- 
able us to reduce the class under consideration to a simpler form. We also present the equivalence 
transformations of the reduced class and a special subclass. 

Note 1. The value of the arbitrary element m is undefined if /i = 0. In this case, it is most 
convenient for one to assume m equal to n + 1. Sometimes the value m = 1 is also acceptable. 

The usual equivalence group G~ of class is formed by the nondegenerate point transfor- 
mations in the space of {t,x,u, f,g,h,n,m), which are projectible on the space of {t,x,u), i.e. 
they have the form 

{i,x,u) = {T\T-,T^){t,x,u), 

if,g,h,h,m) = {Tf,T9,T\T",T"'){t,x,u,f,g,h,n,m), 

and transform any equation from class (^Q) for the function u = u{t, x) with the arbitrary 
elements {f,g,h,n,m) to an equation from the same class for the function u = u{t,x) with the 
new arbitrary elements (/, g, h, n, rh). 

Theorem 1. G~ consists of the transformations 

i = 5it + 82-, x = ip{x), H = 63U, 

where 6j {j = 0,3) are arbitrary constants, SqSiS^, ^ 0, ip is an arbitrary smooth function of x, 
+ 0. 

It appears that class admits other equivalence transformations which do not belong to G~ 
and form, together with usual equivalence transformations, a generalized extended equivalence 
group. Restrictions on transformations can be weakened in two directions. We admit the 
transformations of the variables t, x and u can depend on arbitrary elements (the prefix "gen- 
eralized" [26]), and this dependence are not point necessarily and have to become point with 
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respect to {t,x,u) after fixing values of arbitrary elements. The explicit form of the new arbi- 
trary elements {f,g,h,h,rh) is determined via {t,x,u, f,g,h,n,m) in some non-fixed (possibly, 
nonlocal) way (the prefix "extended"). We construct the complete (in this sense) generalized 
extended equivalence group G"" of class (0) , using the direct method [22] . 

Theorem 2. The generalized extended equivalence group G~ of class (0) is formed by the 
transformations 



t = 5it + 52, X = (f{x), u = ^{x)u, 

/ = ] — nr/) 9 = :o in g, h= ■ — ■ — --rn, n = n, m = m, 

where 6j {j = 0,1,2) are arbitrary constants, 5o5i ^0, <p is an arbitrary smooth function of x, 
ipx ^ 0. The function ip{x) is determined by the formula 

. f(l - (n+ l)F(j;))"^, , ^, . , f dx 

^^^) = |ei^W, n = -l' '^^''^=^'i W)^ '- 

Note 2. The above representation of the function ^/^(x) guarantees that this family of transfor- 
mations is continuously parameterized with the parameter n. 

From theorem 1 we deduce that the transformation 

/dx 
— -, u = u (4) 
9[.x) 

maps equation ^ to 

f{i)ui = (u"ni)£ + /i(x)u'", 

where /(x) = g{x)f{x), g{x) = 1 and h{x) = g{x)h{x). (Likewise any equation of form can 
be reduced to the same form with f{x) = 1.) That is why, without loss of generality, we can 
restrict ourselves to investigation of the equations of the general form 

f{x)ut = (n"u,.)x + h{x)u"'. (5) 

All results on symmetries, solutions and conservation laws of class can be extended to class (pQ) 
with transformations (0J). 

It appears that we can deduce the equivalence group for class © from theorem 1 and 2 by 
setting g = g = 1. The results are summarized in the following theorem. 

Theorem 3. The generalized equivalence group of the class where n / —1, consists of 
the transformations 

t = dit + d2, x=- —=:ip(x), u = djipx^^+'^u, 

~ 3n+4 ~ , , -J m4-3n + 3 

f = didj'^ipx 2"+2/, h = 67''^^'^^ (fx 2n+2 h, h = n, m = m, 

where 5j (j = 1, 7) are arbitrary constants, 6161 7^ and 636^ — 8465 = ±1. 
For n = —1 transformations from group take the form 

i = 6it + 62, X = 63X + 64, u = 5^e^^^u, 

f= — — /j h = — 7^ — ; — h, h = n, m = m, 

where 5j (j = 1,6) are arbitrary constants, 616365 / 0. 
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These equivalence transformations can be employed to simplify the forms of /(x) and h{x) 
in the study of Lie symmetries of © which takes place in the next section. 

Note 3. Theoreminican be reformulated in the following way: The usual equivalence group G~ 
of class (jSJ for values n 7^ — 1 is formed by the transformations 

1 1 

1 1 1 

n = (1 - (n + 1){63X + 64)) "+in = e 2n+2 2^+2 =; ^[x)u, 

~ n 3n4-4 ~ — m-f-JT-4-1 m-|-3n4-3 

/ = 5ie2"+2(^^ 2n+2j^ h = e 2n+2 y,^. 2n+2 /^^ h = U, 171 = 771, 

where 5j (j = 1,4) and e are arbitrary constants, eSiS^ ^ 0. 

This representation of the transformations from is more clumsily than in the theorem 
but it shows that the family of the above transformations is continuously parameterized with 
the parameter n, including the value n = —1. In fact, 

lim <f{x) = - + hm ij{x) = e'53^+^4_ 

n^ — 1 £ £03 n-^ — l 

Since the parameters n and m are invariants of all the above equivalence transformations, 
class (or class can be presented as the union of disjoint subclasses where each from the 
subclasses corresponds to fixed values of n and m. This representation allow us to give other 
interpretation of the above results. For example, the generalized equivalence group from 
theorem ini more exactly its projection to the complementary set of variables and arbitrary ele- 
ments, can be considered as a family of usual equivalence groups of the subclasses parameterized 
with n and m. 

The question occurs: could one obtain wider equivalence groups for some of the subclasses, 
apriori assuming the parameters n and m satisfy a condition? The answer is positive in the case 
of m = n + 1. Namely, the following statement is true. 

Theorem 4. The class of equations 

f{x)ut = {g{x)u^u,)^ + h{x)u^+^ (6) 
admits the equivalence group C^^^^i consisting of the transformations: 
i=5it + 52, x = (p{x), u = ■ip{x)u, 

-J, 9= ,o„i9 g> h = 6Q pr—^ , (n = n), 



where ip and are arbitrary functions of x, 6j {j = 0, 1, 2) are arbitrary constants, 6o6iipxip ^ 0. 

It should be emphasized that is the usual equivalence group of class © even through 

n is assumed as an arbitrary element. Moreover, it is wider than the generalized extended equiv- 
alence group of whole class (0). We will also called G'^^J^_^_^ as the conditional equivalence 
group of class (pQ) under the condition m = n -1- 1 on the arbitrary elements. It is a nontrivial 
conditional equivalence group of class (0) in the sense that it is not a subgroup of G^ . Let us 
note that the notion of conditional equivalence group was first used in classification of systems 
of two non-linear Laplace equations [36] (see also [37] ) . 

Gauging g with the condition g = 1, we impose the restrictions 5o(px = V'^"^^ on parameters 
of G^^^_^_l and obtain the equivalence group G5"^=„+i of the class of equations 

f{x)ut = {u^Ux)x + h{x)u^+\ (7) 
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Corollary 1. The equivalence group Gi^^j^_^i of class Q is formed by the transformations 
i = 6it + 62, X = ip{x), u = ip{x)u, 

f = ^i^f^ ^ = ^0 ' in = ri), 

where if and ip cltc arbitrary functions satisfying the condition JqVx = V'^"'''^; (j = 0) 1> 2) are 
arbitrary constants, doditp 7^ 0. 

Similarly to , Gi^^,^^i is the generalized equivalence group for the whole class © since 
transformations of variables explicitly depend on the arbitrary element n. And be- 
comes the usual equivalence group under restriction to any subclass with a fixed value of n. 

Simultaneous usage of both the generalized equivalence group and the conditional equiva- 
lence group Gi^^^_^i via gauging arbitrary elements leads to crucial simplification of solving the 
group classification problem and presentation of the obtained results. See note Inland example 
for some explanations. 

There exists another condition on arbitrary elements, which gives a non-trivial conditional 
equivalence group. All such conditions and the corresponding conditional equivalence groups 
are systematically investigated as components of classification of admissible transformations in 
section 13 

Note 4. Due to physical sense of equation (pQ), the function u should satisfy the condition u >0. 
In this case we have to demand for the multipliers of u to be positive in all transformations. If 
we avoid positiveness of u then we have to use the modular of u as base of powers which are 
not determined for negative values of base. The same statement is true for similar expressions 
in transformations and other places. The necessary changes in formulas are obvious. 



3 Lie symmetries 

In the section we present a complete classification of Lie symmetries for class © . We search for 
operators of the form T = rit, x, u)dt + i{t, x, u)dx + r]{t, x, u)du, which generate one-parameter 
groups of point symmetry transformations of equations from class ©. These operators satisfy 
the necessary and sufficient criterion of infinitesimal invariance, i.e. action of the r-th pro- 
longation r^*"^ of r to the (r-th order) DE results in identically zero, modulo the DE under 
consideration. Here we require that 

r^^^f{x)ut - u^'uxx - nu'^-^ul - h{x)u'^] = (8) 

identically, modulo equation ©. 

After elimination of ut due to ©, equation (jS)) becomes an identity in five variables, t, x, 
u, Ux and Uxx- In fact, equation (jHJ is a multivariable polynomial in the variables Ux and Uxx- 
The coefficients of the different powers of these variables must be zero, giving the determining 
equations on the coefficients r, ^ and r/. Since equation has a specific form (it is a quasi-linear 
evolution equation, the right hand side of © is a polynomial in the pure derivatives of u with 
respect to x etc), the forms of the coefficients can be simplified. That is, r = T(t), £, = ^{t, x) [22] 
and, moreover, rj = C,{t,x)u and = ^(x) (see Sectional). To derive the latter condition, we 
partially split with respect to u. Finally we obtain the classifying equations which include both 
the residuary uncertainties in coefficients of the operator and the arbitrary elements of the class 
under consideration: 

fxi = f{nQ + Tt-2ix), 2{n+l)Qx=ixx, 

uf\t - u^^^fC. + u'^ihfxi - fhxi + (1 - m)fhC - fhrt) = 0. 
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The third equation is a polynomial in u and should be split with respect to u for all values of the 
parameters n and m, which correspond to the nonlinear case n 7^ 0. Splitting gives essentially 
different results in three exclusive cases: 1. m 7^ 1; 2. m = 1; 3. m = n+1. The obtained 

equations enable us to derive the forms of r(t), ^(x), Citjx), f{x) and h{x) depending on values 
of n and m and consequently the desired Lie symmetries will be constructed. 

In table 1 we list tuples of parameter-functions f{x) and h{x), the constant parameters m 
and n and bases of the invariance algebras in all possible inequivalent cases of Lie symmetry 
extension. The operators from table 1 form bases of the maximal Lie invariance algebras iff 
the corresponding values of the parameters are inequivalent to ones with more abundant Lie 
invariance algebras. In cases 2 and 5 we do not try to use equivalence transformations as much 
as possible since otherwise a number of similar simplified cases would be derived, see remarks 
after the table. 

Note 5. It should be emphasized that we adduce only the cases of extensions of maximal 
Lie invariance algebra, which are inequivalent with respect to Gi if m 7^ n + 1 and with 
respect to Gl^^^_^_l if m = n + 1 or /i = 0. During solving the group classification problem, 
we derive some cases of extensions, which are transformed with equivalence transformations 
to representatives from the table in a non-trivial way. Below we give some examples of such 
transformations . 

Example 1. The 'variable coefficient' equation 

|x| ut = {u'^Ux)x + £\x\ "+1 n™, n 7^ — 1, my^n-\-l, (9) 

is transformed to the "constant coefficient" equation 

% = (10) 

(case 3 of table 1 if m 7^ 1 and case 6 if m = 1) by the transformation 
1 i_ 1 

t = t, X= — =:Lp{x), U = \x\ " + 1 U = 2n+2t( 
X 

which belongs to the equivalence group Gi of the class ©. Hence, using the maximal Lie 
invariance algebra of case 3, we can derive the basis elements of the maximal Lie invariance 
algebra of ©: 

dt, in + l)x^dx + xudu, 2(1 — m)tdt — (n + 1 — m)xdx + ^^—^-^—^—^udu 

n + 1 

if m 7^ 1 and 

dt, e~^"*(9j + Xudu), {n + l)x'^dx + xudu, n{n + l)xdx — {n + 2)udu 

if m = 1. In the latter case we additionally assume n 7^ —4/3 for the algebra to be really 
maximal. 

Example 2. In an analogous way, the 'variable coefficient' equation 

e'ut = {—] +ee'u (11) 

\ U J X 

is reduced by the transformation t = t, x = x, u = e^u from the equivalence group G^ to the 
equation (|10)) with n = —1 and m = 1 (case 6 with n = —1). The basis elements of the maximal 
Lie invariance algebra of Ullj) have the form 

dt, e^\dt + eudu), dx - udu, xdx - {x + 2)udu. 
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Table 1. Results of group classification of class Q 





n 




h{x) 


Basis of A™"" 


General case 


1 


V 


V 


V 


a* 


2 


V 


/l(^) 


fti(a;) 


9t, (d + 2fo — pn)tdt + ((n + l)aa;^ + 63; + c)dx + (aa: + p)udu 


3 


V 


1 


e 


dt, dx, 2(1 - m)t9t + (1 + n - m)xdx + 2udu 


m = 1, ft / 0, (ft//), = 


4 


V 


V 


£/ 




5 


V 


/l(^) 


£f 


dt,e~^"\dt + eudu),n{{n + l)ax^ + bx + c)dx + {nax + 2b + d)udu 


6 


7^ 

' 3 


1 


e 


dt,dx, e"^"'(9t + £udu),nxdx + 2udu 


7 


4 
3 


1 


e 


dt, dx, ei^^{dt + eudv), —^xdx + 2udu, —^x^d^ + xudu 


m = n + l or h — 


8 


V 


V 


V 


dt,ntdt — udu 


9 


^ 3 


1 




dt,ntdt — udu, 2tdt + xdx 


10 


^ 3 


1 


e 


dt,ntdt — udu,dx 


11 


' 3 


1 





dt,dx,ntdt — udu, 2tdt + xdx 


12 


4 
3 




a 


dt,tdt + judu,dx ~ \udu 


13 


4 
3 


1 





dt,dx, \ tdt + udu, 2tdt + xdx, — \x'^dx + xudu 



Here a is arbitrary constant, e = ±1, 

-(3n + 4)aa; + d 



/i(a::) = exp 



(71 + l)ax^ + + c 



dx 



hi (x) — exp 



-(3(n + 1) + in)ax + (1 + m - n)p - 2b 
(n + l)ax'^ + 6a; + c 



dx 



and it can be assumed up to equivalence with respect to that the parameter tuple (a, b, c, d) takes only the 
following non- equivalent values: 

(£,0,1,0) if n = -l or (1,0, l,d') if n / -1, (0, l,0,d'), (0,0,1,1), 

where d' is arbitrary constant. In all the cases we put g{x) — 1. 

In case 8 the parameter-functions / and ft can be additionally gauged with equivalence transformations from 
Gi[m=„4-i. For example, we can put / = 1 if n 7^ —4/3 and / = otherwise. 
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Note 6. Right choice of a gauge of the parameter tuple (/, g, h) is another crucial point of our 
investigations. Really, the major choice have been made from the very outset of classification 
when the parameter- function g was put equal to 1. It is the gauge that leads to maximal 
simplification of both the whole solving and the final results. Although all gauges are equivalent 
from an abstract point of view, only the gauge g = 1 allows one to exhaustively solve the problem 
of group classification of class with reasonable quantity of calculations. Even after the other 
simplest gauge / = 1 chosen, calculations become too cumbersome and sophisticated. 

The expression of transformations from the conditional equivalence group G'^^^_^_l, i.e. the 
equivalence group of the subclass © separated from class (0) with the condition m = n + 1, 
contains one more arbitrary function of x in comparison with the equivalence group of the 
whole class. This fact makes an additional gauge of the parameter tuple (/, g, h) is possible and 
necessary in case of m = n + 1. Different ways of the additional gauge with the conditional 
equivalence group (or are discussed in the next example. 

Example 3. Let m = n + 1 and the gauge g = 1 he already fixed. The kernel A^^^^^^j^i 
of maximal Lie invariance algebras of the subclass is generated by the operators dt and 
ntdt — udu (case 8 of table 1). 

It seems on the face of it that the optimal choice for the additional gauge is /i = 0. This 
gauge results to essential simplification of both the general form of initial equations and the 
determining equations. Indeed, the system of determining equations are then reduced to the 
following one: 

e/x = « - 2^, + r^)/, = 2(n + l)Cx, = 0, Ci = 0. 

Hence, Q = Cix + C2, ^ = Ci(n + + C3X + C4, r = C^t + Cg. The main difficulty arises 
after taking into account the first equation. It implies that extension of is possible iff 

where a, /3, 7 and 5 are arbitrary constants, (a,/3, 7) / (0,0,0). Moreover, the equivalence 
group of the gauged class '5 = 1, /i = 0' is formed by the transformations from which 
additionally satisfy the conditions d^ipx = ip'^'^'^'^, {ip~^"'~^^'^ipx)x = 0. Although classification 
may be carried out in this way, complexity of both the expression of / in some extension cases 
and equivalence transformations to be applied results in necessity of careful and sophisticated 
investigations of reducibility between different cases. 

Another way is to begin with gauging the parameter-function /. The obtained gauge 
4 

/ = 1 if n 7^ and / = or (/, /i) = (1,0) otherwise 

3 

seems to be more complicated than h = 0. But it is the gauge that straight leads to exhaustive 
classification, and formulation of results having a simplest form (cases 8-13 of table 1). 

Note 7. We tested different ways of classification. Implemented 'experiments' result in the 
following conclusion. Application of equivalence transformations is more effective on earlier 
stages of classification. An expedient way is to gauge arbitrary elements as much as possible on 
preparatory classification stage and additionally to gauge arbitrary elements under classification 
every time when sufficient conditions to do it arises. The worst choice is to refuse to use gauges 
completely. 

After analyzing the classification adduced in table 1, we can formulate the following theorem 
which is similar to a one from [38]. 

Theorem 5. If an equation of form (OJ) is invariant with respect to a Lie algebra of dimension 
not less than four then it can be reduced by means of point transformations to a one with constant 
values of the parameters f, g and h. If m ^ l,n + l, the similar statement is also true for three- 
dimensional Lie invariance algebras. 



9 



4 Additional equivalence transformations 



Let a class of differential equations be classified with respect to its equivalence group G~ and pro- 
jection of G~ on the space of the equation variables be narrower than the whole pseudo-group 
of point transformations in this space. Then there can exist point transformations between 
C^^-inequivalent cases of symmetry extension, which are called additional equivalence trans- 
formations. Knowledge of them is important since it simplifies further application of group 
classification results. 

The first non-trivial and interesting example of such transformations was given in [24] (see 
also [14, Chapter 10]). Under classification of 'constant coefficient' nonlinear reaction-diffusion 
equations Q in [9,10], equations of the form 



ut 



where b is an arbitrary constant, arose as cases of symmetry extension, which are inequivalent 
with respect to the corresponding equivalence group for different values of b. It was shown in 
[24] that the transformation 



= 'I X — \ U — \ 

[ tan(Y^-6/3x), 6 < ' [ cos^{y^ -b/3x)u, b<0 

maps any such equation with the equation 

which is of the same form with b' = 0. Therefore, it is an additional equivalence transformation 
in class (jSJ. It is an additional equivalence transformation in classes and © and, at the 
same time, a member of the equivalence group Gl^^^_^_l of class ((TJ. That is why the case 6/0 
does not arise under our classification. Let us note that application of generalized, extended and 
conditional equivalence groups simultaneously with the usual one results in implicit but effective 
usage of additional equivalences. 

Another additional equivalence transformation in class (jSJ [14, Chapter 10]) 

t' = —e'^\ x' = x, u' = e-''u (12) 
en 

links the equations 

Ut = {u^Ux)x + su and u[, = {u''^u'^,)^,. (13) 

This transformation belongs to no equivalence group found in section |2I It can be obviously 
extended to a wider subclass of Namely, it connects the equations 

f{x)ut = {g{x)u"'Ux)x + £f{x)u and f{x')u[, = {g{x')u'^u'x>)x' 

and, therefore, reduces cases 4-7 of table 1 to the set of cases = n + 1 or h = 0', i.e. it is an 
additional equivalence transformation in class Any other additional equivalence transfor- 
mations is a composition of transformations of form (|12)1 and transformations from GJ"^^^.,.]^. 
(It is proved in the next section in the framework of admissible transformations.) For example, 
if n = — I equations H13() are also connected by the mapping 

t' = — e^"*, x' = -, u' = e-^^x^u, 5 = ±1. 
en X 
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Seeming scarcity of additional equivalence transformations in class e.g. in comparison 
with the class f{x)ut = [g{x)A{u)ux]x + h{x)B{u)Ux admitting a large number of additional 
equivalence transformations [16,38], is connected with usage of conditional equivalence groups 
under classification. As it is demonstrated in the next section, the transformations of form H12|l 
can be also included in the framework of conditional equivalence. 

As a result, the following statement is derived. 

Theorem 6. Up to point transformations, a complete list of extensions of the maximal Lie 
invariance group of equations from class is exhausted by cases 1-3 and 8-13 of table 1. 



5 Classification of form-preserving (admissible) transformations 

Due to special structure of equations from class the following problem can be solved com- 
pletely. To describe all point transformations each of which connects a pair of equations from 
class (^. Such transformations are called form-preserving [22] or admissible [37] transforma- 
tions. See [37] for stronger definitions. They can be naturally interpreted in terms of the category 
theory [41]. Let us note that there exists an infinitesimal equivalent of this notion [6]. 

Since class can be gauged with transformations from its usual equivalence group it is 
enough for one to solve the similar problem in class ©. To do it, we consider a pair of equations 
from the class under consideration, i.e. © and 

= {u^ux)x + h{xW (14) 

and assume that these equations are connected via a point transformation of the general form 

t = T{t,x,u), x = X{t,x,u), u = U{t,x,u). 

We have to derive determining equations in functions T, X and U and to solve them depending 
on values of arbitrary elements in and (|14j) . 

After substitution of expressions for the tilde- variables into ()14() . we obtain an equation in the 
tildeless variables. It should be an identity on the manifold determined by © in the second-order 
jet space over the variables {t, x\u). Splitting of this identity with respect to the unconstrained 
variables in a true way results at first in the equations 

— — Xu — — 

that agrees with results on more general classes of evolution equations [22,38,41]. We split 
further, taking into account the above equations and also splitting with respect to u partially. 
As a results, the following relations are deduced: 

h = n, fh = m or (m, m) G {(1, n + 1), (n + 1, 1)}, 

T = T{t), X = X{x), U = V{t,x)u, TtXxV^O 

2(n + l)XxVx = XxxV, 

fXx' 



yn 



f Tt 



f ^ i...m . I^t^_ 1 {^rn^x \ n+1 



^-/i«™ + /^n = — (^F'^^J U^+' + hV^U^. (15) 

Splitting of equation ()15() with respect to u and subsequent integration of the determining 
equations depends on values of m and rh appreciably. Consider different cases separately. 
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If 771 = m 7^ 1, equation (|15j) implies = hence Tu = 0, i.e. V = V{x) and T = Sit + 62- 
The other conditions are 

/ h yn+i-m 

m = m/l,n+l: '7 = T 



f 'X,' h X,2 ' 
n = -l: X^^ = 0, (InF)^.^. = 0, 



Therefore, in the case m ^ l,n + 1 (the case m = n + 1) any admissible transformation is 
determined by a transformation from the equivalence group Gi (the conditional equivalence 
group 

If m = m = 1, the function V is conveniently presented as y = \Tt\~^^""il^{x). The value 
n = — 1 is again singular: 

n = -1 : = 0, (In = ^ X = dgx + ^4, = ^se^^"^, 

n/-l: (IX.I-V2) 0, 2(n + l)^ = ^ ^ X = ^5:^, ^ = 

where (^siJe — 6465 has a fixed value which can be assumed equal to ±1. The arbitrary elements / 
and h are transformed by the formulas 

f ip"" h 1 h 1 f 1\ ^ fh\ 1 fh 

■ sign Tt, — = — — H — TfT 5 hence 



In view of the latter condition, {h/ f)x / and Tu = if {h/ f)x / 0. Therefore, any admissible 
transformation in this case corresponds to a transformation from the equivalence group . 

Different situation is in case of {h/ f)x = 0. Then {h/f)x = 0, i.e. the condition {h/f)x = is 
invariant with respect to G^. Let the constants h/f and h/f be denoted by a and a, respectively. 
Then we get the following equation for the function T{t): 

^\ 1 ~ 
— = —na h na. 

TtJt Tt 

We integrate these equation and present the general solution in such form that continuous 
dependence of it on the parameters a and a is obvious: 

^naT _ ^nat _ e"°"^ — 1 

aa^O: = Si \- S2, a = 0, a^O: = Sit + S2, 

na na na 

^nat _ ^ 

a / 0, a = 0: T = Si h S2, a = a = 0: T = Sit + S2- 

na 

The above-described transformations form the equivalence group G^^^^ (/i//)^=o sub- 
class of class which is separated with the conditions g' = 1, m = 1, {h/f)x = 0. It is a 
generalized equivalence group even through n is fixed since it contains transformations with 
respect to t, which depends on (constant) ratio of the arbitrary elements h and /. Unlike 
Gi ,^^^^i, we do not use G'^^^^ ^^^-^ under group classification of class © since applica- 
tion of this conditional equivalence group does not have crucial influence on classification and 
the corresponding conditions on arbitrary elements are less obvious. At the same time, trans- 
formations from G^^^j^ (h/f)x=o P^^^ '■^^^ °^ additional equivalence transformations after 
completing the classification (see the previous section). 
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If /i/i 7^ 0, m = 1 and m = n + 1 then the determining equations yield 

h iTu . /h\ ^ , 1 



a; 



Hence both equations © and (|14j) are mapped to the subclass '/i = 0' of the class under 
consideration. These mappings are realized by the transformations from the corresponding 
conditional equivalence groups. We can assume that their images coincides. Therefore, the 
admissible transformation is a composition of a transformation 7i of © to the equation f{x)ut = 
{vP'Ux)x and transformation T2 of the equation f{x)ut = {u"'Ux)x to (|14|) . The transformation 7i 
belongs to G^^^^ ^^^^^ and the transformation T2 belongs to GJ^^^^^^, and we can put 
X = X, ijj = 1 in Ti. 

The case hh ^ 0, m = n + 1 and m = 1 is considered in similar way. 

All the possible cases are exhausted. 

Let us summarize investigation of the admissible transformations in class © in the following 
statement. 

Theorem 7. Let the equations f{x)ut = {u^Ux)x + h{x)u'^ and f{x)ut = {u^Ux)x + h{x)u"^ be 
connected via a point transformation in the variables t, x and u. Then 

n = n and either m = m or (m, m) = (1, n + 1) or (m, m) = (n + 1, 1). 

The transformation is determined by a transformation from 

a) Gi if either m 7^ 1, n + 1 or m = 1, {h/f)x 7^ 0; 

b) G5;;^=„+i ifm = m = n + l; 

c) G~^^-^ ifm = m=l, {h/f)x = 0; then also {h/f )x = 0. 

If m = 1 and m = n + 1 then {h/f)x = and the transformation is composition of two 
transformations from G^^^-^ (h/f) =0 ^^'^ ^im=n+i ^^^^ intermediate equation having h = 0. 
The case m = n + 1 and m = 1 is similar. 

Corollary 2. Class © with n ^ can be presented as the union of normalized subclasses 
separated by the conditions 

h^O, m^l,n + l; m = 1, {h/f)x^O; m = 1, {h/f)x = 0; m = n + l. 

Only two latter subclasses have a non-empty intersection, and the intersection being the normal- 
ized subclass 'h = 0'. 

Roughly speaking, the class of differential equations is called normalized if any admissible 
transformation in this class belongs to its equivalence group. See [37] for strong definitions. 



6 Basic definitions and statements on conservation laws 

In this section we give basic definitions and statements on conservation laws, following the well- 
known textbook of Olver [30] in general outlines. Then we formulate the notion of equivalence 
of conservation laws with respect to equivalence groups, which was first introduced in [39]. This 
notion is a base for modification of the direct method of construction of conservation laws, which 
is applied in section [7] for exhaustive classification of local conservation laws of equations . 

Let £ be a system L{x, U(^p^) = of I differential equations = 0, . . . , = for m unknown 
functions u = (n^, . . . , u™) of n independent variables x = {xi, . . . , Here U(p) denotes the set 
of all the derivatives of the functions u with respect to x of order no greater than p, including 
u as the derivatives of the zero order. Let denote the set of all algebraically independent 
differential consequences that have, as differential equations, orders no greater than k. We 
identify £(fc) with the manifold determined by in the jet space 



13 



Definition 1. A conserved vector of the system C is an n-tuple F = {F^{x, . . . , F^{x, U(^r))) 
for which the divergence DivF := DiF^ vanishes for all solutions of C (i.e. DivF|^= 0). 

In Definition n and below Di = D^^ denotes the operator of total differentiation with respect 
to the variable Xj, i.e. Di = 5^. + u'^^du^, where u'^ and u'^^ stand for the variables in jet 
spaces, which correspond to derivatives d^°'^u°'/dxi^ . . .dx"" and du'^/dxi, a = («!,...,«„)) 

G N U {0}, \a\: = ai + ■ ■ ■ + On- We use the summation convention for repeated indices and 
assume any function as its zero-order derivative. The notation means that values of V are 
considered only on solutions of the system C 

Definition 2. A conserved vector F is called trivial if F^ = F^ + F^, i = l,n, where and 
are, likewise F*, functions of x and derivatives of u (i.e. differential functions), F* vanish on the 
solutions of £ and the n-tuple F = {r, ■ ■ ■ , F"^) is a null divergence (i.e. its divergence vanishes 
identically) . 

The triviality concerning the vanishing conserved vectors on solutions of the system can 
be easily eliminated by confining on the manifold of the system, taking into account all its 
necessary differential consequences. A characterization of all null divergences is given by the 
following lemma (see e.g. [30]). 

Lemma 1. The n-tuple F = {F^, . . . , F"), n >2, is a null divergence (Div F = 0) iff there exist 
smooth functions v^^ (i,j = l,n) of x and derivatives ofu, such that v^^ = —v^^ and F^ = Djv^^ . 

The functions u'-' are called potentials corresponding to the null divergence F. If n = 1 any 
null divergence is constant. 

Definition 3. Two conserved vectors F and F' are called equivalent if the vector- function F' — F 
is a trivial conserved vector. 

The above definitions of triviality and equivalence of conserved vectors are natural in view 
of the usual "empiric" definition of conservation laws of a system of differential equations as 
divergences of its conserved vectors, i.e. divergence expressions which vanish for all solutions 
of this system. For example, equivalent conserved vectors correspond to the same conservation 
law. It allows us to formulate the definition of conservation law in a rigorous style (see e.g. [43]). 
Namely, for any system C of differential equations the set CV(£) of conserved vectors of its 
conservation laws is a linear space, and the subset CVo(/2) of trivial conserved vectors is a 
linear subspace in CV(£). The factor space CL(£) = CV(£)/ C\q[C) coincides with the set of 
equivalence classes of CV(£) with respect to the equivalence relation adduced in Definition 

Definition 4. The elements of CL(>C) are called conservation laws of the system £, and the 
whole factor space CL(£) is called as the space of conservation laws of C. 

That is why we assume description of the set of conservation laws as finding CL(C) that is 
equivalent to construction of either a basis if dimCL(£) < oo or a system of generatrices in 
the infinite dimensional case. The elements of CY{C) which belong to the same equivalence 
class giving a conservation law are considered all as conserved vectors of this conservation 
law, and we will additionally identify elements from CL(>C) with their representatives in CV(£). 
For F G CV(£) and £ CL(£) the notation F € will denote that F is a conserved vector 
corresponding to the conservation law J^. In contrast to the order rp of a conserved vector F as 
the maximal order of derivatives explicitly appearing in F, the order of the conservation law T 
is called Yi\m.\rp \ F G .F}. Under linear dependence of conservation laws we understand linear 
dependence of them as elements of CL(£). Therefore, in the framework of "representative" 
approach conservation laws of a system L are considered linearly dependent if there exists linear 
combination of their representatives, which is a trivial conserved vector. 
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Let the system C be totally nondegenerate [30]. Then application of the Hadamard lemma 
to the definition of conserved vector and integrating by parts imply that divergence of any 
conserved vector of C can be always presented, up to the equivalence relation of conserved 
vectors, as a linear combination of left side of independent equations from C with coefficients A'^ 
being functions on a suitable jet space 

DivF = A''L''. (16) 

Here the order k is determined by C and the allowable order of conservation laws, ^ = 1, L 

Definition 5. Formula ((TB|) and the /-tuple A = (A^, . . . , A') are called the characteristic form 
and the characteristic of the conservation law Div F = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of C Since C is nondegenerate, 
the characteristics A and A satisfy (jlGj) for the same F and, therefore, are called equivalent iff 
A — A is a trivial characteristic. Similarly to conserved vectors, the set Ch(£) of characteristics 
corresponding to conservation laws of the system £ is a linear space, and the subset Cho(£) of 
trivial characteristics is a linear subspace in Ch(£). The factor space Chf(£) = Ch(i2)/ Cho(£) 
coincides with the set of equivalence classes of Ch(£) with respect to the above characteristic 
equivalence relation. 

We can essentially simplify and order classification of conservation laws, taking into account 
additionally symmetry transformations of a system or equivalence transformations of a whole 
class of systems. Such problem is similar to one of group classification of differential equations. 

Proposition 1. Any point transformation g maps a class of equations in the conserved form 
into itself. More exactly, the transformation g: x = Xg{x,u), u = Ug{x,u) prolonged to the 
jet space J^^^ transforms the equation DiF^ = to the equation DiFg = 0. The transformed 
conserved vector Fg is determined by the formula 

D ■ X- ■ 1 

Fg(x,-Lt(r)) = ^^'^1 F^(x,U(^)), i.e. Fg{x,U(r-)) = ^^^-^{D^x)F{x,Uf^r)) (17) 

in the matrix notions. Here \Dxx\ is the determinant of the matrix D^x = (D^jXi). 

Note 8. In the case of one dependent variable (m = 1) g can be a contact transformation: 
X = Xg{x,U(^i)), tt(i) = iig(i)(x, n(i)). Similar notes are also true for the statements below. 

Definition 6. Let G be a symmetry group of the system C Two conservation laws with the 
conserved vectors F and F' are called G -equivalent if there exists a transformation g ^ G such 
that the conserved vectors Fg and F' are equivalent in the sense of Definition |31 

Any transformation g £ G induces a linear one-to-one mapping g^, in CV(£), transforms 
trivial conserved vectors only to trivial ones (i.e. CVo(-C) is invariant with respect to g^^) and 
therefore induces a linear one-to-one mapping g{ in CL(i2). It is obvious that gf preserves linear 
(in)dependence of elements in CL(£) and maps a basis (a set of generatrices) of CL(£) in a 
basis (a set of generatrices) of the same space. In such way we can consider the G-equivalence 
relation of conservation laws as well-determined on CL(£) and use it to classify conservation 
laws. 

Proposition 2. // the system C admits a one-parameter group of transformations then the 
infinitesimal generator X = + ri^du<^ of this group can be used for construction of new 
conservation laws from known ones. Namely, differentiating equation jj?} ) with respect to the 
parameter e and taking the value e = 0, we obtain the new conserved vector 

r = -X^r)F' + (DjOF^ - {Dji^)F\ (18) 

Here X(^) denotes the r-th prolongation [30, 32] of the operator X . 
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Note 9. Formula p8j) can be directly extended to generalized symmetry operators (see, for 
example, [4]). A similar statement for generalized symmetry operators in evolutionary form 
= 0) was known earlier [13, 30]. It was used in [20] to introduce a notion of basis of 
conservation laws as a set which generates a whole set of conservation laws with action of 
generalized symmetry operators and operation of linear combination. 

Proposition 3. Any point transformation g between systems C and C induces a linear one- 
to-one mapping g^ from CV(£) into CV(£), which maps CYq{C) into CVo(>C) and generates a 
linear one-to-one mapping gf from CL(£) into CL(£). 

Corollary 3. Any point transformation g between systems C and C induces a linear one-to-one 
mapping gi from Chf(£) into Chf(£). 

It is possible to obtain an explicit formula for correspondence between characteristics of C 
and C. Let C^' = A'^'"/:^, where A^'^ = A^^^^L*", A'^'"" are differential functions, a = (ai, . . . , a„) 
runs the multi-indices set (oj EN U {0}), /i, = 1, 1. Then 

\^' = M'^'*{\D^x\~y). 

Here A'^'^* = (— Z))" • K^^^^ is the adjoint to the operator A^'*. For a number of cases, e.g. if C 
and C are single partial differential equations (/ = 1), the operators are simply differential 
functions (i.e. A'^'^" = for |a| > 0) and, therefore, A'^^* = A'^'^. 

Consider the class C\s of systems Cq: L{x,Ui^p-^,9{x,U(^p^)) = parameterized with the 
parameter-functions 9 = 9{x, U(^p-j). Here L is a tuple of fixed functions of x, u^p-^ and 9. 9 denotes 
the tuple of arbitrary (parametric) functions 9{x, u^p-^) = {9^{x, U{p)), • • • , 9^{x, Ui^p-^)) running the 
set S of solutions of the system S{x,U(^p-j,9(^q-^{x,U(^p-^)) = 0. This system consists of differential 
equations on 9, where x and ii(p) play the role of independent variables and 0^^) stands for the 
set of all the partial derivatives of 9 of order no greater than q. In what follows we call the 
functions 9 arbitrary elements. Denote the point transformations group preserving the form of 
the systems from £[5 as G~ = G~(L, S). 

Consider the set P = P{L,S) of all pairs each of which consists of a system Cq from £[5 
and a conservation law of this system. In view of Proposition |21 action of transformations 
from G~ on C\s and {CV{Co) \ 9 £ S} together with the pure equivalence relation of conserved 
vectors naturally generates an equivalence relation on P. 

Definition 7. Let 9,9' e S, T € CL{Ce), G CL(£0/), F ^ T , F' ^ P. The pairs {Ce,^) 
and {Cqi,J-') are called -equivalent if there exists a transformation g G G~ which transforms 
the system Cq to the system Cqi and such that the conserved vectors Fg and F' are equivalent 
in the sense of Definition |21 

Classification of conservation laws with respect to C will be understood as classification 
in P with respect to the above equivalence relation. This problem can be investigated in the 
way that is similar to group classification in classes of systems of differential equations, especially 
it is formulated in terms of characteristics. Namely, we construct firstly the conservation laws 
that are defined for all values of the arbitrary elements. (The corresponding conserved vectors 
may depend on the arbitrary elements.) Then we classify, with respect to the equivalence group, 
arbitrary elements for each of that the system admits additional conservation laws. 

In an analogues way we also can introduce equivalence relations on P, which are generated 
by either generalizations of usual equivalence groups or all admissible point or contact transfor- 
mations (called also form-preserving ones [22]) in pairs of equations from C\s- 

Note 10. It can be easy shown that all the above equivalences are indeed equivalence relations, 
i. e. they have the usual reflexive, symmetric and transitive properties. 
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7 Conservation laws 



We search (local) conservation laws of equations from class applying the modification of 
the most direct method, which was proposed in [39] and is described in the previous section. It 
is based on using the definition of conservation laws, the notion of equivalence of conservation 
laws with respect to a transformation group and classification up to the equivalence group of 
a class of differential equations. In view of results of the previous section, it is sufficient for 
exhaustive investigation if we classify conservation laws of equations only from class © • Let us 
note that other kinds of the direct methods, which are based e.g. on the characteristic form of 
conservation laws [1,2] could be also applied. 

There are two independent variables t and x in equations under consideration, which play a 
part of the time and space variables correspondingly. Therefore, the general form of constructed 
conservation laws will be 

DtF{t, X, U^r)) + D^G{t, X, U(^r)) = 0, (19) 

where Dt and Dx are the operators of total differentiation with respect to t and x. The com- 
ponents F and G of the conserved vector (F, G) are called the density and the flux of the 
conservation law. Two conserved vectors {F,G) and {F',G') are equivalent if there exist such 
functions F, G and H oit, x and derivatives of u that F and G vanish for all solutions of C and 
F' = F + F + DxH, G' = G + G- DtH. 

The following lemma on order of local conservation laws for more general class of second-order 
evolution equations, which covers class Q, is used below. 

Lemma 2 ( [18]). Any local conservation law of any second- order (1 + 1)- dimensional quasi- 
linear evolution equation has the first order and, moreover, there exists its conserved vector with 
the density depending at most on t, x, and u and the flux depending at most on t, x, u and Ux- 

Theorem 8. A complete list of equations |3|) having nontrivial conservation laws is exhausted 
by the following ones 

1. m = n + 1. 

n/-l: {ip'fu, -(p'u''ux + ipi^), (p\ i = l,2. 

n = — 1: {xfu, —xu~^Ux-\-lnu), x; {fu, —u~^Ux),l- 

2. m = 1, h = i^if. 

n/-l: {xe-f'^fu, e-''\-xu''ux + ^)), xe-f"*; (e"^*/?^, -e"''*u"M^), e"/^*. 

n = — 1: {xe~^^fu, e~'^^{—xu^^Ux + lnu)), xe"'^*; (e~^*/M, —e~^^*u~'^Ux), e~^^. 

Here f3i and /?2 are arbitrary constants. The functions = ^^{x), i = 1,2, form a fundamental 
set of solutions of the second-order linear ordinary differential equation ifxx + (?t- + l)hip = 0. 
(Together with constraints on the parameter-functions f and h we also adduce conserved vectors 
and characteristics of the basis elements of the corresponding space of conservation laws.) 

Proof. In view of lemma 121 we can assume at once that F = F{t,x,u) and G = G{t, x,u,Ux). 
Let us substitute the expression for ut deduced from © into H19() and decompose the obtained 
equation with respect to Uxx- The coefficient of Uxx gives the equation u'^f~^Fu + = 0, 
therefore G = —u'^f~^FuUx + G{t,x,u). Taking into account the latter expression for G and 
splitting the rest of equation (|19|) with respect to the powers of Ux, we obtain the system of 
PDEs on the functions F and G of the form 

F„„ = 0, -n"f^) +G„ = 0, Ft + -u^Fu + Gx = ^. (20) 



17 



Solving first two equations of ^ yields F = x)fu + x), G = Ju"^ + x) 
The case n = — 1 is special with respect to integration of the function n". So 



G 



-^u^^Ux + ^x^^u + G^{t,x), n = -l. 



(It is convenient to separate / as a multiplier in the coefficient of u in the expression of F.) In 
further consideration the major role is played by a differential consequence of system (|2()|) that 
can be written as 

f^t + + mh<^u"'-^ = 0. (21) 

Indeed, it is the unique classifying condition for this problem. In all the classification cases we 
obtain the equation Fj^ + = 0. Therefore, we can assume -F'^ = G'^ = up to conserved 
vector equivalence and have to suppose for existence of non-trivial conserved vectors that ^ ^ 0. 

It follows from (|21l) that equation © possesses non-trivial conserved vectors only for the 
special values of the parameters n and m, namely, if 

(n, m) G {(0, 0), (0, 1), (n', n' + 1), (n', 1), n' G M}. 

We exclude from consideration the cases (0, 0) and (0, 1) since they correspond to the linear case 
of equation Q . 

1. m = n + 1. Equation H21|) is split to two equations $t = and ^xx + {n + l)h^ = 0, i.e. 
$ depends only on x and runs the set of solutions of the second equation which is homogeneous 
linear second-order ordinary differential one. ^ = Cix + C2 if n = — 1. 

2. m = 1. Split of ()21() results in the equations ^xx = and f^t + = 0. The equa- 
tion ^xx = gives <^ = <I>^(t)x + <I>^(t). Substitution ^ in this form to the second equation yields 
^Ifx + -|- ^^hx -|- ^^h = 0. It follows from the latter equation in view of the condition 
$ / that /i = /i/ and $ = G3e~^*x + Qe"^*. □ 

Corollary 4. For any (nonlinear) equation from class dimension of the space of local 
conservation laws equals to either or 2. In the second case the equation can he reduced by 
a point transformation to an equation from the same class, which has the values g = I and 
h = 0. Then a basis of the corresponding space of characteristics is formed by the identically 
unity (ip^ = Ij o-nd the function (p^ = x. 



8 Similarity solutions 

The Lie symmetry operators found as a result of solving the group classification problem can 
be applied to construction of exact solutions of the corresponding equations. The method 
of reductions with respect to subalgebras of Lie invariance algebras is well-known and quite 
algorithmic to use in most cases; we refer to the standard textbooks on the subject [30,32]. 
Exhaustive group analysis, including Lie reductions, of nonlinear constant-coefficient reaction- 
diffusion equations with nonlinearities of general form was carried out in [10] (see also [14]). 
A number of exact solutions of such equations are tabulated (see e.g. [35]). Then the problem 
on Lie reductions and Lie exact solutions of equations from class , which have four- or five- 
dimensional Lie invariance algebra, can be assumed as solved in view of theorem [21 

The Lie reduction algorithm is more effective for equations having nice symmetry properties. 
Together with the above facts, it gives solid argumentation for us to choose the cases from table 1 
with at least one non-constant parameter-function of x and a three-dimensional Lie invariance 
algebra. Namely, these are cases 9 and 12. (Let us remember that case 5 can be transformed 
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by an additional equivalence transformation to the case with the same value of / and e = and 
then reduced by equivalence transformations to one of cases 9-12 depending on the values of n 
and /.) 

We have shown that equations 

Ut = {u'^Ux)x + ax"^!*""*"^ and e^ut = {u~3Ux)x + (xu~3 

(table 1, case 9 and case 12) admit the three-dimensional Lie invariance algebras generated by 
the operators 

Xi = dt, X2 = tdt du, X3 = xdx H 9„ 

n n 

and 

3 3 

Xi = dt, X2 = tdt + Xs = dx- -udu 

correspondingly. Due to strong simplification of these cases with equivalence transformations, 
the operators are simplified to simple linear combinations of translation and scale operators and, 
therefore, are handy for further usage. The both tuples of operators satisfy the commutation 
relations 

[Xi , X2] = Xi , [Xi , X3] = 0, [X2 , X3] = 0. 

It means that the algebras are isomorphic to the algebra 02 © 0i being the direct sum of the 
two-dimensional non-Abelian Lie algebra 02 and the one-dimensional Lie algebra 0i. An op- 
timal set of subalgebras of 02 0i can be easily constructed with application of a standard 
technique [30,32]. Another way is to take the set from [33]. (In this paper optimal sets of sub- 
algebras are listed for all three- and four-dimensional algebras.) The used optimal set consists of 

one-dimensional subalgebras: (X2 — /UX3), (^3), {Xi^X^), (Xi); 
two-dimensional subalgebras: {Xi , X3 — 1^X2) , {Xi , X2) , 

where /x and ly arc arbitrary constants. 

Lie reduction to algebraic equations can be made only with the first two-dimensional sub- 
algebra; the second one does not satisfy the transversality condition [30]. The corresponding 
ansatzes and reduced algebraic equations have the form. 

9.u = Cx'', where a = (7"+^((n + 1)^^ - a + a) = 0; 

12. M = Ce^^, where (7 = --(i/ + l); C^^ (^^ - 3a) = 0; 

Here C is an unknown constant to be found. The reduced equations are compatible and have 
non-trivial (non-zero) solutions only for some values of a and, moreover, become identities for 
these values of a. As a result, the following stationary solutions are constructed: 

9. u = Cx'^ , where (n -|- 1)(7^ — a + a = 0; 
12. u = Ce""", where = 3a; 

Here C is an arbitrary constant. These solutions can be also obtained with step-by-step reduc- 
tions with respect to one-dimensional subalgebras. 

The ansatzes and reduced equations corresponding to the one-dimensional subalgebras from 
the optimal system are collected in table 2. 
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Table 2. Similarity reductions of cases 9 and 12 



N 


X 




u = 


Reduced ODE 


9.1 


X2 - ^iX3 




t 1 ^(y) 




9.2 


X'i 


t 


x~ lp{lo) 


nVc^ ~ (aw^ + 2n + 4)v3"+^ = 


9.3 


X-i±X^ 






iv^fu,)^ ± i.Ji/'^ ^1(^4- auj-^ip"+^ = 


9.4 


Xi 


X 


Lp{uj) 




12.1 


X2 - flXi 


X + ptXrit 




{(p~^ifiu^)u, - ^J.e'^fu^ - |(/i + l)e'^(y5 + aV'~* =0 


12.2 


X-i 


t 




le^Po; + (3 - 16q)<^'4 = 


12.3 




x^^t 






12.4 


x^ 


X 




. 4 . 1 

(V3 s'/'")" + 3=0 



Some reduced equations are integrated completely: 

1 

4\ 



9.2. (p 



9.4. 



C - (an + 2 + -juj 



if a' = 0, 

if a' > 0, 
1 1 

uj 2("+i) ^Ci sin((T In w) + C2 cos((t lna;))"+i, if a' < 0, 



Ciw2("+i) (Inw + C2)"+i, 



where a' = 1 — 4a(n + 1), >ci^2 

12.2. : 



1 ± 



a' 



Cl(a; + C2)"^ if « = 0, 

12.4. = <{ (Cie'^^ + C2e-^'^)"^ if a > 0, 

(Ci sin(cra;) + C2 cos(crw)) ^, if a < 0, 
where x = a/3, cj = -y/— a/3; 



The above solutions of reduced equations results after substitution to the corresponding ansatzes 
to exact solutions of the initial equations. 

Let us note that the obtained solutions can by transformed by equivalence or admissible 
transformations to exact solutions of much more complicated reaction-diffusion equations than 
the investigated ones. 

9 Conclusion 

In this paper different properties and objects concerning variable coefficient (l+l)-dimensional 
nonlinear reaction-diffusion equations are investigated in the framework of classical group 
analysis: equivalence transformations, Lie symmetries, conservation laws and exact solutions. 

The cornerstone of presented investigation is application of different generalizations of the 
usual equivalence group, including generalized, extended and conditional equivalence groups 
and admissible transformations, in an optimal way via gauging arbitrary elements. It is the tool 
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that makes principal simplifications in solving the group classification problem, construction of 
conservation laws and finding exact solutions. The set of admissible transformations of class 
has an interesting structure which is exhaustively described by theorem [3 

It seems natural for one to consider other classes of PDEs in the same way. The direct 
extension is investigation of variable coefficient (l+l)-dimensional nonlinear reaction-diffusion 
equations for different (e.g. non-power) nonlinearities. Thus, we will complete studying the 
case n = which is omitted in this paper. 

The adduced results form a basis for advanced analysis of class with 'modern' symmetry 
methods. In particular, classification of conservation laws creates the necessary prerequisites for 
studying potential symmetries [5]. It is undoubted that generalized equivalences will continue 
to play a singular role in these investigations. 
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